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11,6  design  and  selection  of  appropriate  fender  systems  must  consider  the  dynamic  I 
interaction  between  the  ship  and  port  or  pier  structure.  Energy  absorption 
characteristics  of  marine  fender  systems  vary  as  a  function  of  fender  geometry 
material,  load  time  history,  load  spatial  distribution,  local  hull  stiffness  and 
frequency  of  loading.  The  dynamic  interaction  of  marine  fender  systems  with 
ship’s  hull  is  studied  in  time  domain.  The  frequency  dependence  of  the 
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1.  INTRODUCTION 


In  order  to  achieve  more  appropriately  designed  marine  fender  systems,  an 
understanding  of  the  vessel-fender  dynamic  interaction  is  essential.  The 
dynamic  analysis  can  describe  more  accurately  the  fender's  energy  absorption 
characteristics  and  operational  performance  requirements. 

There  are  several  approaches  which  consider  the  total  energy  of  a  berthing 
ship  to  be  absorbed  by  the  fender/pier  system  in  addition  to  the  magnitude 
of  the  fender  reaction  force  generated.  A  simple  and  commonly  used  approach 
is  the  energy  method;  e.g.,  Lee  [10],  Brolsma  et.al.  [1],  Accurate  predic¬ 
tions  using  this  method  involve  knowledge  of  hydrodynamic  coefficients  and 
system  stiffnesses  for  the  total  system.  Empirical  values  are  commonly  used 
in  -this  method .  System  damping  has  been  ignored  in  this  approach.  Alternate 
statistical  approaches  consider  information  obtained  from  model  measurements 
to  determine  the  design  energy  value  of  the  berthing  impact  and  fender  energy 
absorption;  e.g.,  Sevendsen  [13].  With. the  design  risk  selected,  the  design 
value  of  the  energy  to  be  absorbed  by  the  fender  system  can  be  determined. 

The  disadvantage  of  this  method  is  in  obtaining  the  necessary  statistical 
information  for  problem  solution.  Kim  [?]  proposes  an  approximate  and 
simpler  method.  The  idea  is  based  on  the  use  of  time  average  of  the  kinetic 
energy  of  the  berthing  ship  during  the  time  interval  of  fender  compression. 
This  method  is  applicable  only  to  the  ship  in  calm  water  and  berthing 
broadside. 


The  time  domain  solutions  of  forces  and  motions  have  been  developed  by 
van  Oortmerssen  [15]  and  Fontijn  [3].  Both  of  these  methods  use  Impulse 
Response  Function  techniques.  It  is  easy  to  adopt  different  external 
forces^ and  some  other  factors  in  the  time  domain  method.  The  procedure 
presented  herein  is  based  upon  a  simplified  convolution  integral  for  added 
mass  and  damping  calculations,  which  eliminates  the  disadvantage  of  the 
.expensive  calculation  in  the  above  time  domain  solutions.  Although  the 
mathematical  problem  is  formulated  herein  in  sway,  yaw,  surge  and  roll 
motions,  only  the  lateral  motion  results  are  presented  here.  The  hydro- 
dynamic  coefficients j  as  functions  of  frequency,  can  be  determined  theo¬ 
retically  using,  two-dimensional  strip  theory  or  three-dimensional  source 
distribution  method.  In  order  to  verify  the  current  technique  with 
Fontijn's  results,  the  same  hydrodynamic  coefficients  are  used  herein. 


T  ** 


2.  DYNAMIC  BESPONSF  OF  THF  SHIP  AND  THE  BEPTHING  5TPUCTUPF 


During  the  berthing,  the  ship  will  undergo  dynamic  motion  which  has  six 
degrees  of  freedom;  namely,  swaying,  yawing,  rolling,  surge,  heave  and 
pitch.  The  heave  and  pitch  motions  are  of  little  consequence  in  energy 
dissipation  and  may  be  neglected.  In  the  following  analysis,  it  is 
assumed  that  there  is  no  sliding  contact  along  the  fender’s  surface. 

Als p,  the -coupling  effects  between  each  mode  have  been  neglected. 

Consider  the  dynamic  equilibrium  of  the  center  of  gravity  of  the  ship  as 
shown  in  Figure  1,  the  equations  of  motion  are: 


(mk  *  ak)  *k  +  bk*k  "  fk  k  "  1’2,4,6 

where  x1f  x2,  xj,  and  xg  are  the  surge  motion  (x*  positive  forward), 

•  sway  notion  (y,  positive  to  port),  roll  -motion  ($),  and  yaw  motion  (9), 
respectively.  The  mk,  ak  and  bk  are  inertia  mass,  hydrodynamic  mass 
and  hydrodynamic  damping  in  the  corresponding  directions.  fk  represents 
-the  external  fender  reaction  force  on  the  ship  in  the  k  direction. 

For  the  berthing  structure  with  the  effective  mass  Ms,  the  dynamic  equations 

are.: 


Vs  +  Vs  +  Vs 


-f, 


(2) 


My  +  y  y  +  ■  -f  0 

s  s  y  s  y  s  2 


(3) 


where  K  and  vx  y  are  the  structure  stiffness  and  damping,  respectively. 
The  displacement  of  the  point  of  contact  can  be  calculated  after  the  dynamic 
equations  (1)'  to  (3)  are  solved. 
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,} j£j  If  the  elastic  deformation  of  the  ship  hull  is  under  consideration,  let  Kh, 

'  !?j  and  S^  and  y  jj  define  the  stiffness,  deflection  and  the  damping  coeffi¬ 

cients  of  the  hull  at  the  point  of  contact  respectively.  The  equations  of 

P:  motion  of  the  ship  hull  are: 
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“h^hx  +  vhx\x  +  KhxShx  *  fl 

Vhy  +  Vhy  +  Vhy  m  f2 


(4) 


(5) 


in  the  case  of  a  very  rigid  berth,  the  deflection  of  berthing  structure  can 
be  neglected  due  to  the  high  stiffness  of  the  structure.  The  effects  of  the 
ship’s  local  stiffness  can  ;be  investigated  with  this  method. 
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The  solution  of  dynamic  equations  is  carried  out  by  the  numerical  integration 
method  with  appropriate  choice  of  problem  parameters. 


3i  HTDKODYHAKEC  MASS  AND  DAMPING 

The  hydrodynamic  mass  and  damping  are  important  parameters  to  be  considered 
in  the  determination  of  the  berthing  forces.  The  hydrodynamic  msss  is 
governed  by  the  following  factors: 

-  Ship  characteristics  (beam,  draft,  size) 

Under-keel  clearance 

—  Water  depth 

-  Type  of  berthing  structures  (open,  semi-open,  solid) 

-  Fender  characteristics 

-  Berthing  modes 

-  Berthing  velocity. 
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Most  of  the  methods  to  determine  the  hydrodynamic  characteristics  do  not 
take  into  account  all  of  the  above-mentioned  factors.  A  recent  review  by 
Kray  [9]  has  presented  a  good  summary  of  the  state-of-the-art  on  hydro- 
dynamic  mass  determination.  The  most  common  practice  is  to  have  a  constant 
added  mass  as  presented  by  Vasco  Costa  [16],  Komatsu  and  Salman  [8]. 

Hayashi  and  Shirai  [1*3  presented  a  theoretical  formulation  of  ship  added  mass 
in  shallow  water  as  a  function  of  the  ratio  of  draft  to  water  depth,  the 
Froude  number  of  the  ship  and  the  coefficient  of  head  loss  of  the  counter 
flow  under  the  hull.  Oda  [12]  conducted  experiments  to  verify  their  theory. 

In  order  to  have  better  agreement  in  very  shallow  water  berthing,  Oda 
modified  the  theory  by  neglecting  the  shear  stress  acting  on  the  hull  and 
also  used  a  friction  parameter  as  a  function  of  sway  force  coefficient,  con¬ 
traction  coefficient  and  water  depth  to  draft  ratio.  The  theory  modified  by 
Oda  will  be  used  to  compare  with  the  existing  experiments  and  present  method. 

The  constant  added  mass  coefficient  is  not  appropriate  for  berthing  of  ships 
in  shallow  water  with  different  berthing  speeds.  Oortmerssen  [13]  and  Endo  [2] 
calculated  the  hydrodynamic  characteristics  of  ship  as  function  of  frequencies 
in  ihe  shallow  water  based  on  three-dimensional  approach.  Comparisons  have 
been  made  with  the  results  for  cases  available  in  the  literature.  Generally 
speaking,  the  agreement  is  good.  In  the  study  by  Fontijn  [3],  the  two- 
dimensional  hydrodynamic  characteristics  at  low  frequencies  were  modified 
to  fit  the  experimental  data.  Time  histories  of  fender  force  in  Fontijn’s 
paper  compared  well  with  the  present  method. 

4.  EFFECTS  OF  LOCAL  HULL  STIFFNESS 

In  the  design  and  selection  of  the  most  appropriate  fender  system,  considera¬ 
tion  should  not  only  take  into  account  the  energy  absorption  capacity  of  the 
fender  systems  but  also  the  strength  requirement  of  ship  side  structure  due 
to  berthing  impact.  The  bending  rigidity  of  hull  structure  in  the  horizontal 
direction  is  extremely  high,  so  that  the  deformation  from  the  whole  ship, 
assumed  as  a  beam,  is  not  critical.  Two  kinds  of  deformation  are  considered 
by  Kawakami  et.al.  [5,6].  There  are  the  deformations  of  hull  plate  between 
frames  and  the  deformations  of  hull  plate  between  frames  and  trnasverse  bulk¬ 
heads.  The  analytical  treatment  of  the  hull  deformations  and  stresses  due 
to  berthing  impact  is  also  presented  in  the  report  edited  by  Matsunami 
et.al.  [11],  In  their  analysis,  the  deflections  and  stresses  of  ship  hull 
are  considered  as  a  static  reaction. 


UM 


When  the  contacting  region  between  ship  hull  and  fender  is  small,  the  load 
distribution  within  the  contacting  area  can  be  generally  considered  to  be 
uniform.  In  the  present  dynamic  analysis,  the  constant  contact  area  on  the 
ship  hull  is  either  on  the  center  of  the  plate  or  on  the  beam,  but  can  be 
simulated  anywhere  once  described. 


5.  MATHEMATICAL  METHOD 


As  mentioned  earlier,  a  basic  difficulty  arises  in  the  time  domain  solution 
dealing  with  the  frequency-dependent  hydrodynamic  coefficients.  The  force 
derived  from  drag,  and  added  inertia  is  a  complex  function  of  ship  movement 
at  each  given  frequency.  The  time  history  of  a  force  denoted  by  f(t)  can  be 
described  by  means  of  the  convolution  integral  as: 


K(t--r)  y,(r)  dr 


Where  y(x)  is  the  time  history  of  the  ship  motion,  e.g.,  sway.  The  Kernal 
function  K(t)  is  the  Fourier  transformation  of  the  complex  transfer  function 
of  F(u}  with  respect  to  the  ship  motion,  which  is  given  by: 


w 

/ 


F(w)  e*Ut  dw 


This  procedure  has  been  used  in  a  number  of  investigations  to  obtain  the 
histories  of  forces,  e.g.,  Fontijn  C3l  and  van  Oortmerssen  [15]. 


In  the. present  application,  the  convolution  operation  is  implemented  by 
digital  computation  in  a  different  manner.  Tou  [ 14]  has  shown  that,  for  a 
single  input  sinusoidal  frequency  of  amplitude  a,  the  instantaneous  value 
of  the  output  time  sequence,  f(nAt),  can  be  predicted  by  a  weighted  sequence 
of  n  previous  time  history  steps  as  follows: 


—  «  7  «  •  *  *  «  "v  .  .  v  *  •  -  * 

t.-  '.-  .-  v 


n 

f(nAt)  ■  a  £  wm  sin  [w  (n-o)  At]] 
in  “  1 


(8) 


where: 

wn  s  weighting  coefficients 
At  =  time  interval 

However,  the  instantaneous  values  of  this  function  can  also  be  given  in 
terms  of  an  input  sinusoid  as: 


f(nAt)  *  a  sin  (wnAt  +  a)  ’(9) 


Where  a  is  the  phase  angle  of  the  output  relative  to  the  input.  Equations 
(8)  and  (9)  set  equal  to  each  other  give  a  single  equation  relating  the 
transfer  function  F  to  the  weighting  coefficients  wm  for  a  given  fre¬ 
quency  w  .  When  dealing  with  inputs  containing  a  number  of  frequencies, 
this  method  leads  to  a  set  of  simultaneous  equations  for  the  value  of  wm 
in  terms  of  the  values  of  the  transfer  function  at  each  frequency.  These 
equations  can  be  solved  by  using  a  least  squares  approach  for  a  set  of 
weighting  coefficients,  wQ.  When  applying  this  method,  the  operation 
is  over  the  past  time  history  of  the  ship  motion.  The  convolution  integral 
is  then  replaced  by  a  truncated  summation: 


n 

f(nAt)  «  53  wm  F  £(m-n)Atj| 

m  «  1 


(10) 


After  the  evaluation  of  the  added  inertia  force  and  drag  force,  the  dif¬ 
ferential  equations  are  solved  by  Runge^Kutta  step-by-step  integration 
procedure. 


6.  RESULTS 


The  analytical  technique  described  herein  has  been  applied  to  the  case  in 
Fontijn's  paper  [3].  The  main  dimensions  of  the  ship  used  in  the  model  are: 
waterline  length  =  2.^438  m,  beam  =  0.375  m,  draft  =  0.15  m,  block  coefficient  = 
1.0,  ship  mass  s  137. 24  kg,  water  depth  =  0.2  m.  The  added  mass  and  damping 
coefficients  were  determined  by  Fontijn  using  2-D  theory  with  three-dimensional 
effects  on  the  low  frequency  range.  The  same  hydrodynamic  coefficients  are 
used  in  order  to  have  a  fair  comparison.  Two  types  of  fenders  are  considered 
in  the  centric  berthing:  a  linear  fender  of  constant  stiffness,  k^,  equal 
to  1 4*0  kg/m,  and  a  nonlinear  fender  force  represented  by 
k1  Ayf  +  k2(  Ayf-df),  in  which  k^  =64  kg/m,  kg  =  113  kg/m,  df=0.00664  m 
and  Ay^  is  the  deflection  of  fender. 

Figure  2  shows  the  comparison  between  present  method  with  Fontijn's  calculated 
and  measured  fender  forces  as  a  function  of  time.  This  analytical  method  gives 
good  agreement  with  the  experiment.  Hydrodynamic  coefficients  have  great  influ¬ 
ence  on  the  fender  response.  The  use  of  the  existing  results  is  the  first  step 
in  the  development  of  this  method.  The  nonlinear  fender  stiffness  is  calcu¬ 
lated  for  the  same  ship  in.  Figure  3.  Fender  force  and  ship  motions  for 
eccentric  berthing  on  linear  fenders  are  given  in  Figures  4  and  5.  The  ship 
characteristics  and  water  depth  are  same  as  centric  impact.  Fender  stiffness, 
klf  equals  to  65  kg/m.  In  these  figures,  the  quantity  mQ  is  defined  as 
p  LBDmg(  pLEDe02+  mg)"*1.  Theoretically,  the  program  has  been  formulated 
to  accept  characterization  of  any  fender  system  for  the  calculation 
of  energy  absorption.  Various  generic  fender  system  algorithms  have  been 
developed  which  will  be  part  of  further  program  development. 

Hayashi  and  Shirai’s  method  [4]  is  considered  a  simple  way  to  calculate 
the3 added  mass  coefficient.  They  took  into  account  the  shallow  water,  ship 
speed  and  characteristics.  That  method  was  modified  by  Oda  [12]  to  have  a 
better  agreement  with  experiments.  Based  on  the  modified  theory,  fender 
force  on  the  same  model  used  in  this  paper  is  presented  in  Figure  6.  The 
approaching  phase  shows  better  agreement  than  the  detaching  phase;  Oda 
shows  the  similar  behavior  in  his  results. 

<If  one  takes  into  consideration  the  ship  hull,  then  the  fender  force  depends 
oh  the  hull  stiffness  and  mass.  For  the  model  considered  above,  we  assume 
the  local  hull  is  simulated  by  aluminum  plate  with  size  0.3  &  by  0.3  m  and 
0.006  »  thickness;  total  weight  on  the  plate  is  15  kg.  The  fender  reaction 
is  almost  never  affected  by  the  hull  under  these  circumstances  (stiffness  is 
about  100,000  kg/m).  If  the  hull  stiffness  is  reduced  to  1000  kg/m  (academic 
purpose).,  then  the  maximum  fender  force  changes  about  6.5  percent. 


7.  CONCLUSIONS 


The  frequency  dependence  of  the  hydrodynamic  coefficients  can  be  easily 
represented  by  a  simplified  convolution  integral.  This  method  has  good 
agreement  with  existing  theory  and  experiment.  The  hydrodynamic  coeffi¬ 
cients  significantly  affect  the  vessel-fender  interaction.  More  experi¬ 
mental  and  theoretical  studies  are  needed  in  order  to  include  additional 
environmental  effects.  Local  hull  effects  do  not  significantly  affect 
the  dynamic  response  problem.  This  method  provides  a  means  of  analyzing 
the  dynamic  ship/fender  problem  in  a  simplified  cost-effective  manner. 


8.  ACKNOWLEDGEMENTS 

This  work  has  been  supported  by  the  Office  of  Naval  Research,  Contract  No. 
NOOOl!|-82-C-OOl8,  under  the  direction  of  Dr.  Eugene  A.  Silva.  The  authors 
also  wish  to  express  their  appreciation  to  Drs.  T.  T.  Lee  and  J.  G. 
Giannotti  for  their  many  valuable  comments  and  continuous  encouragement. 


u 

•  J 


r 


j*  C.V. >.  J*- » i?  •  -  •  ». 


!-**<•■&?• «?.W  »'*  -.  %  „L  -"  _  w* '.  *■  *  V*  ?.'■%  •„*,  V».T  «  M 


'•  B 
'  id 


m 


m  n 

;:>•’  at.’ 


Jj  ?3j 
»•  Si 


•V 

9;«g 

>4  ft 
IV  ft! 


IS 
1  W 


9.  REFERENCES 

1.  Brolsma,  J.  V.,  Hirs,  J.  A.,  et.al., "On  the  Fender  Design 
and  Berthing  Velocities, "  24th  International  Navigation  Congress, 

1977. 

2.  Endo,  H.,  "Motions  of  a  Vessel  Moored  Hear  Vertical  Walls  and 
with  Limited  Under-Keel  Clearance,"  Thesis,  Univeristy  of  Hawaii, 

1981. 

3.  Fontijn,  H.  L.,  "The  Berthing  of  a  Ship  to  a  Jetty,"  Journal  of  the 
Waterway  Port  Coastal  and  Ocean  Division.  ASCE,  Vol.  106,  No.  WW2, 

May  1980. 

4.  Hayashl,  T.  and  Shirai,  M.,  "Added  Masses  of  Large  Tankers  Berthing 
to  Dolphins,"  ASCE,  15th  Coastal  Engineering  Conference,  Honolulu, 
Hawaii,  Vol.  Ill,  1976. 

5.  Kawakami,  M.,  et.al.,  "On  the  Relations  of  Fender  to  Strength  of  Ship 
Side  Structure,"  Hiroshima  University,  Japan  Research  Report,  Vol.  24, 
Nc.  i;  1975. 

5.  kawakami,  M.,  et.al.,  ?0n  the  Required  Function  of  Fender  Considering 
of  Prevention  of  Damage  of  Ship  Side  Structure,"  Hiroshima  University, 
Japan  Research  Report,  Vol.  24,  No.  2,  1976. 

7.  Kim,  C.  H.,  "Sway  Added  Mass  and  Impact  Load  of  a  Large  Tank  to  a 
Jetty  in  Shallow  Water,"  The  Arctic-Offshore-Deepsea  Systems  Sympo¬ 
sium,  ASME,  1982. 

8.  Komast3U,  S.  and  Salman,  A.  H.,  ^Dynamic  Response  of  the  Ship  and  the 
Berthing  Fender  System  After  Impact,"  Proceedings  of  JSC-E,  No.  222, 
April  1972. 

.9.  Kray,  C.  J.,  "Safety  of  Ships  and  Structures  During  Berthing  and 
Mooring,"  Journal  of  Marine  Technology  Society,  Vol.  16,  No.  1, 

1982. 

10.  Lee,  T.  T.,  "A  Design  Criteria  Recommended  for  Marine  Fender  Systems," 
proceedings  of  the  11th  International  Conference  on  Coastal  Engi¬ 
neering,  London,  ASCE,  1968. 


•j 


11.  Matsunami,  H.,  et.al.,' "Design  of  Fender  Systems.  Working  Group  oh 
Fender  System  Design,*  Japanese  National  Section  of  PIANC,  March  1980. 

12.  Oda,  K.,  "Virtual  Masses  of  Supertankers  Berthing  to  Dolphins," 
Proceedings  of  the  28th  Japanese  Conference  on  Coastal  Engineering, 
1981. 

13.  Svendsen,  I.  A.,  "The  Statistical  Approach,"  PIANC  International 
Study  Commission  for  Improving  Fender  Design  Systems,  Edinburgh, 

1981. 

14.  Tou,  J.  T.,  Digital  and  Sampled-Data  Control  Systems.  McGraw-Hill 
Book  Company,  1959. 

15.  van  Oortmerssen,  G.,  "The  Motions  of  a  Moored  Ship  in  Waves,"  Thesis, 
the  Delft  University  of  Technology,  The  Netherlands,  1976. 

16.  Vasco  Costa,  F.,  "The  Berthing  Ship,  the  Effect  of  Impact  on  the 
Design  of  Fenders  and  Other  Structures,"  The  Dock  and  Harbor 
Authority,  1964. 


0  1  2 


4  5  6  7 


t  *  k^/m 


TIME  HISTORY  ol  FENDER  FORCE:  LINEAR  FENDER 


FIGURE  2 


TIME  HISTORY  of  FENDER  FORCE: 
NON-LINEAR  FENDER 


FIGURE  4 

TIME  HISTORIES  or  ECCENTRIC  BERTHING 
eo/L  =  0.167 


nOn 


v/^h-.0293 


4fc 


>  1.2 


FIGURE  5 

TIME  HISTORIES  ol  ECCENTRIC  BERTHING 
ep/L=  0.5 


>  8 

N 

E° 

%  " 

•P  A 


.*0 

r'o^ - 


o 

o  _ 


I . 


o 

E  1.4 


_  Theory 

O  Exp  ... 

>Fon1ijn 
.  Theory! 

v/^gh*.0256 


fel 


\0 


JBi 


k7/rn0 


■O  O  O.*  Ol  «■ 


TIME  HISTORY  of  FENDER  FORCE: 

I  - - 


COMPARISON 


_  Present  Method 

.  Hayashi  et  al 


